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SUMMARY 


O O 
| are 
The electron distribution in the —C—O—P group of the carboxyl-phosphate molecule 
Oo— 
has been calculated by a semi-empirical method. Six electrons are taken into consideration 
explicitly, viz. two z-electrons of the C=O double bond, the lone pair at the bridge O-atom and 
two electrons of the phosphoryl group, corresponding to the usual z-approximation. For com- 
parison, the electron distribution in the carboxyl ion, in the phosphate ion and in polyphosphates 
of the same type as adenosine diphosphate and triphosphate have been calculated. It appears 
from the calculations that the charges due to the “z-electrons” on the C- and P-atoms respec- 
tively are approximately the same for all the molecules considered. The charges on the bridge 
O-atoms in carboxyl-phosphate and in the polyphosphates are on the other hand very different 
from the charges on the other O-atoms of the same group and also from the charges on the O-atoms 
of the carboxyl and phosphate ions. The possibility of explaining some experimental data, by 
means of the results obtained, concerning energy transfer in the living cell is discussed. 


1. Introduction 


Qi. aa 


Carboxyl-phosphate, R—-C—-O—P , can be considered as a representative of 


: OH 

a group of compounds in the living cell which make possible a certain type of energy 
transfer. These compounds are called high-energetic since relatively much energy is 
released when they are hydrolyzed. Kalckar [1] and Oesper [2] have pointed out the 
significance of ‘opposing resonance”’ for their high energy content; e.g. in carboxyl- 
O O 

| Il 7 


phosphate, opposing resonance in the —C—O—P __ group is significant. However, 


: 


as far as the author is aware, no quantum chemical calculation of the electron distri- 
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bution has been performed previously. As such a calculation could be of interest and 
of importance for the understanding of the energy transfer in the living cell, the 
author has made an approximate calculation of the electron density at the different 
O O 
O= 
atoms of the —C—O—P group. For comparison, the electron distribution is 


= 
calculated for the carboxyl ion and phosphate ion as well as for polyphosphates of the 
types, adenosine diphosphate (ADP) and adenosine triphosphate (ATP). These cal- 
culations are reported in the present paper. Some implications concerning the energy 
transfer have been discussed in.two other papers [3, 4], in which preliminary results 
of the present calculations were used. 


2. Method 


The semi-empirical method used here is the same as that used by Fischer [5]. 
Only six of the electrons involved could be taken into consideration explicitly, viz. 
two z-electrons of the C-—O double bond, the lone pair at the bridge oxygen atom 
and two electrons of the phosphoryl group PO;. This approximation corresponds to 
the usual z-approximation. Consequently we will call the dislocalized orbitals of the 
carboxyl-phosphate “‘z-orbitals”, the quotation marks indicating that the orbitals 
are not z-orbitals in the usual sense. 

The calculation has been carried out by the molecular orbital method. Atomic 
orbitals (AO’s) have been adopted for all the atoms of the group and molecular 
orbitals (MO’s) have been approximated by linear combination of atomic orbitals 
(LCAO-MO). An MO ’; over a group of the atoms A, « ... can be expressed as 


¥i= 2 aPas (1) 


where gy, is a normalized AO of the Ath atom and the cj, are parameters. Normaliza- 
tion of the MO’s gives 


favs Vjdt= > cy Cin Sin = (Us (2) 
A. 
where Say = i Pi PudT. (3) 
Then the energy of an electron in the jth MO is 
e= {YH Vaz, (4) 


where H is a one-electron Hamiltonian. The parameters c,, are as usual determined by 
minimizing the energy (4), which gives the condition 


| Ha. —€ Si,|=0 (5) 
where Ay= [oi Hp, dt. (6) 


Each root ¢ of equation (5) corresponds to a molecular orbital, i.e. to a set of c-para- 
meters. 
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Table 1. Ionization energies in eV, taken from Skinner and Pritchard [6]. 


I, is the ionization energy required to remove a p electron from atom A, the atom being in the 
ground state. 


Ig 11.26 Ig-lo —2.35 

Io 13.61 Ip-Ig —2:63 

Ip 10.98 Ig-Ip + 0.28 

We now introduce the following notations 

Hy= OO, iy = Bays (7) 
Vine Bau am 4 Sau (a, at Cay (8) 
01, — hy = 204g (9) 
oo Ez. (10) 


The index C in (10) indicates that the quantity is referred to a carbon atom. Introduc- 
ing (7)—(10) into equation (5) we obtain 


Z2+2 016 Vig t Sas (z a dict Ox¢) 0 
V12 =r Sas (z oe O1c + 02c) z+2 Osc V 23 sl Sea (z 55 O20 on O3c) 0) 
0 Vos + Sos (2 + doc + Jac) 2+ 26s = () 


(11) 


According to the semi-empirical method used here, the values of the integrals 
64, and y,;, should be determined from their connection with experimental data. 
Inserting these values in (11), this equation in z can be solved. The MO’s can then 
be obtained from a calculation of the sets of c-parameters, one set corresponding 
to each z-value. 

The experimentally determined quantity that most closely corresponds to the 
integral 6;, may be half the difference between the ionization energies of the atoms 
and A, but there is no exact correspondence. Upon deliberation, however, it seemed 
that |26,,| should be somewhat smaller than |J,,—J,|. Therefore it appeared neces- 
sary to make several calculations, differing only in the 6,, values used, to obtain an 


Table 2. 26,, values (in eV) used in the calculations. 


Calculations denoted A B C | D 
2860 ~1.0 —1.4 —1.6 —2.0 
: 206p —1.4 -1.8 -1.8 —2.2 
: 9 9 
: 26 + 0.4 + 0.4 + 0.2 + 0.2 
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insight into the dependence of the results on the 6,,, values. Ionization energies given 
by Skinner and Pritchard [6] are listed in Table 1 and the 6,, values used in the 
calculations are listed in Table 2. 

The quantities y,, have been obtained from experimentally determined bond 
energies in the following manner. 

Let us consider a localized MO extending only over the two atoms 4 and yw. The 
distance between the atoms / and yu is assumed to be the same as in the group with 
dislocalized orbitals since it is for this distance that we wish to calculate y,,,. Hence, 
we are considering a hypothetical case. When equation (11) is solved for this case, 
we obtain 


C= Orc Onc + b= ap, Shy ae aC V; Vu tC Bd Ob. ( (1 - a aN: (12) 


Here only the root with the positive sign is of interest since this root corresponds 
to the MO of lowest energy. 

To obtain a quantity that can be calculated from experimental data, it is con- 
venient to form the expression D’ = a, + «, —2¢. Then from (9), (10) and (12) we 
obtain 


D'= to. a 1 sy, Saude i V Vi ae 05, (1 eT S7.,)]- (13) 


D'Sx: feast 
This gives Van = = f+ V2 Ohne (14) 


for the determination of y,,,. Of the two values, the negative one corresponds to the 
definition of y,, in (8). 

To determine y;, from (14), the values of D’, S;, and 6,, must be known. In a 
treatment of z-electrons, D’ is approximately equal to the difference between the 
dissociation energies of a double bond and a single bond, the interatomic distance in 
both cases being the same, viz. the distance to which the y,, and S,, are referred. 
In the calculation of y,, for a o-bond, D’ is approximately equal to thte dissociation 
energy of a o-bond with the interatomic distance to which y;,, and S), are referred. To 
calculate the dissociation energy D, for an arbitrary distance r between the atoms 
when the dissociation energy D, for the equilibrium distance 7) is known, the Morse 
curve 


D,—D,=D, (1-7 °? [15] 


has been used. Here a = Vk/2D,, k being the force constant. Unfortunately all the 
necessary experimental values for the force constants k and dissociation energies 
D, could not be found in the literature. Cottrell’s tables [7] and the values reported 
in Fischer [5] have been used for the estimation of the values used in the present 
calculation. These latter values are listed in Table 3 and the D, values calculated 
from (15) are given in Table 4. 

To obtain values for the overlap integrals S;,, the explicit form of the AO’s, gy, 
must be fixed. The AO’s used in the calculations are given in section 3. When possible, 
the values of the overlap integrals have been obtained by interpolation in Mulliken’s 
tables [8]. The remaining values have been calculated by means of Rosen’s A and 
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Table 3. Values used for the calculation of D,. 


O-C O-P 


Single bond | Double bond | Single bond | Double bond 


ro in A 1.43 1.21 1.66 1.39 
Dissociation energy Dy, Keal/mole 80 156 80 156 


Force constant k, dyne/em 4.4-10° 12.0- 10° 4-10° 11-105 


B functions, which have recently been tabulated by Flodmark [9] for an extensive 
range of parameter values. All of the S;, values obtained are given in Table 5. 

D, values from Table 4, S;, values from Table 5 and 6;,values from Table 2 have 
been inserted into equation (14) for the calculation of y,,. The negative values of 
yi. thus obtained are given in Table 5. 

The y;, values calculated in this way have been inserted into equation (11) together 
with values for 6,, and S),. This equation was then solved by means of the Swedish 
electronic computer Besk, thus giving z-values and the parameters c,,. 

The energy values (4) were obtained from the z-values according to (10). The value 
for %,_ was assumed to be —11.26 eV according to [6], see Table 1. 

The charges on the atoms, arising from the electrons considered in the present 
calculation of MO’s, have been computed by means of Mulliken’s definition [10]. 
For this reason, two populations, which are defined by Mulliken in the following 
manner, have been calculated. The partial gross population NV(j;A) in a singly occupied 
MO, 'V’;, and AO, ,, is defined by 


N (53 Ay =a t+ D> cacy Say. (16) 

7 
The summation over has been carried out over only the nearest neighbours of the 
atom /. The gross atomic population N(A) on the atom A, expressed here for doubly 


occupied space orbitals since there are 2 electrons in each orbital considered, is 
defined by 


N (A)=2>.N (934) =2> Gh +2> > cn ein Sau, (17) 
u J a 


Table 4. Dissociation energies calculated from values in Table 3 according to (15). 


Distances = | Distances 

between C on + 1 K hes le || between F Keal/mole eat as le 

and O, in A pal/mole eal/mo | and O, in A al me {mo 

| 

1.21 — 6 1.45 154 61 
1.23 156 61 1.50 148 70 
1,28 152 70 1.55 142 76 
1.40 136 80 1.60 134 79 
1.43 131 — 
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Table 5. Calculated values for Sj, and yj, where yj, is calculated according to (14) 


and given in eV. 


Values used for the orbital constant ¢, and for dju (in eV) are given. For further information 


about AO’s used, see section 3. 


Distance between C and O, in A Ue} | 1.23 | 1.28 | 1.40 
Soc (2 p%, 2p 7) 0.294 | 0.284 | 0.260 | 0.213 
Yoo(2p%, 2p) for 
2059= = 1.0 — 2.75 — 2.58 =2.17 =—137 
28og= - 1.4 — 2.69 = — 2.10 = 
205q=1—1.6 — 2.65 — — 2.05 ao 
2060= = 2.0 — 2.56 oo —1.94 — 
Distance between P and O, in A 1.45 1.50 1.55 
Sop (2p 9, 3po) for 
pS br = 0.319 0.321 
Ges = 1.60 0.290 — 0.305 
Sop (2p %, 3pm) for 
Cp=1.72 - ou ne 
Cp = 1.60 _ = 0.218 
Sop (2p, 4p 7) for 
Cp = 0.43 — 0.037 0.038 ¢ 
Yop (2p, 3p0) for 2659p = — 1.4 and 
Cp=1.72 = — 1.83 ~2.01 
€p = 1.60 hoe — —1.99 
Yop (2p%, 3p) for 2659p = — 1.4 and 
Cp=1.72 ns = ae 
Cp= 1.60 Bs = — 1.56 
Yop (2p, 4p) for €,=0.43 and 
PMepeet— olky! — — 1.63 —1.30 
2dop= —1.8 —_ —1.52 -1.17 
2dop= — 2.2 = — 1.38 —0.97 


0.038 


where the summation over 7 should be carried out over the occupied space orbitals 
and where the summation over y is as in (16). The charge Q(A) on the atom A is defined 


as 


Q (A) = Ny (A) — N (A), 


(18) 


where N,(A) is the number of electrons in the analogous orbitals of the neutral atom. 
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3. Choice of AO’s and the geometry of the molecules 


In the present application of the method described above to the “‘z-electrons”’ 
of the carboxyl-phosphate, several calculations have been made which differ in the 
orbitals adopted for the phosphoryl group and in the estimate of interatomic 
distances. The different assumptions are collected in Tables 6 and 7 and will now 
be briefly discussed. 

For the O-atoms and the C-atom 2p, hydrogen-like orbitals have been adopted 


S 


5 
mp2 ret? Y (0, 9). (19) 


The ¢-values used are 1.88 for oxygen [11] and 1.56 for carbon [5]. 
For the P-atom, the following kinds of orbitals have been used: 
(a) in calculations I, II, [V and V, a 3p analytical orbital 

ee 


fr? 
ne: ret" ¥ (9, 9), (20) 


where € = 1.60 in calculations I and II and € = 1.72 in IV and V, and 
(6) in calculations III and VI, a 4p analytical orbital 


a 
Pir 3 aeehee Y (8, Q); (21) 
where ¢ = 0.43. 

The ¢-values for the P-atom have been calculated according to Slater. The value 
1.60 refers to the shielding of a 3p electron when the configuration is 1s? 2s72p%3s?3p3, 
the value 1.72 to a 3p electron when the configuration is 1s?2s?2p§3s!3p?3d! and the 
value 0.43 to a 4p electron when the configuration is 1s?2s?2p%3s13p34p!. 

The spatial location of the different atoms is given in Fig. 1. The O,-atom, the 
C-atom, the O,-atom, i.e. the bridge atom, and the P-atom are assumed to be in 
the xy-plane. The atomic orbitals used for these atoms are p, orbitals. The Og and the 
three O, atoms are situated at the corners of a tetrahedron around P. 

In the calculations II, IV and V, it has been assumed that the three O-atoms 
denoted y form a group called (O,),. For this group, the orbital 


2 1 
Por-V a5 p (2 po,1) ae {yp (2, po,2) +p (2 po,s)} (22) 
3 2 


has been used. 2pa,,, represents a 2p orbital of O,,, in the direction of the o-bond of 
P-O,,,. The orbital (22) is of the same symmetry as the orbitals of the atoms C, O, 
and O; when P-O,, is in the yz-plane (see Fig. 1a). In the other calculations, only 
one of the O,-atoms, O,,, has been included in the calculation. In calculation I, O,, 
is assumed to be situated in the yz-plane and, in calculations [III and VI, in the xy- 
plane (see Fig. 1b). The 2p, orbital has been used for O,,. 

The calculations for the carboxyl ion, the phosphate ion and the polyphosphates 
of the types, adenosine diphosphate (ADP) and adenosine triphosphate (ATP), 
have been made using the same parameter values as in certain of the above-mentioned 
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a 


Fig. la Fig. 16 


Fig. 1. Oz, C, Og and P are in the wy-plane. (a) P—O,, is in the yz-plane. Calculations I, II, IV 
and V. (b) P—O,, is in the awy-plane. Calculations III and VI. 


Op aRC, 

Adenosine —P—0,.—P—OH ADP 
on On 
Dies jupte Oath eae 


| I | 
Adenosine —P—O,,—P—O,,—_P—OH ATP 


| | 
OH OH OH : 


Fig. 2. The notations of the atoms in adenosine diphosphate and adenosine triphosphate. 


calculations for carboxyl-phosphate. The parameter values for all the different mole- 
cules are summarized in Table 6. 

The interatomic distances in the calculations have been estimated from known 
distances for pure single and double bonds and from distances in the carboxyl ion 
and the phosphate ion. The distances used in the calculations are listed in Table 7. 


4. Results 


In Table 8, the gross atomic population N(A) according to (17) and the charge 
Q (A) according to (18) due to the “‘z-electrons” for each atom are given. The summa- 
tion over 7 in (17) must be taken over a different number of MO’s for the different 
molecules: (i) for carboxyl-phosphate and adenosine diphosphate, over the three 
lowest orbitals, (11) for the carboxyl ion and the phosphate ion, over the two lowest 
orbitals, and (iii) for adenosine triphosphate, over the five lowest orbitals. In the 
calculation of Q(A), the V(A)-value used is equal to 1 in all cases except for the bridge 
O-atoms, i.e. O; in the carboxyl-phosphate, O,, in ADP and O,, and O,, in ATP. 
For these bridge O-atoms, NV,(A) is equal to 2 since a lone pair, which belongs to the 
neutral atom, is considered. 

Table 9 shows the calculated energy values in eV. The ¢-values were obtained 
from equation (10), the value for %, being assumed to be —11.26 eV. 

Table 10 in the Appendix contains the partial gross population N(j;A) defined in 
(16). 
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Table 6. General summary of atomic orbitals adopted in the calculations. 


Calculations I-VI refer to carboxyl-phosphate. Ox, C, Og and P are situated approximately in 
the ay-plane. For the notations concerning the atoms, see Figs. 1 and 2. Orbital exponent denoted 


by ¢. 
| | Ox Cc | Og | iz | Ov | (O3),, 
Orbital 2p 2p 2p 3p 2p,, o-bond 
pearak c 1.88 | 1.56 | 1.88 | 1.60 | “1.88 = 
Orbital As in cale. I 3 pz Y,,, of (22) 
vp tae c 1.60 a “1.88 
= Orbital ¢ 4p 2 p,, m-bond 
Calc. ITI, VI t = 0.43 Z 1.88 — 
ae Orbital = 3p @,, of (22) 
rd Rl ¢ 1.72 nN 7") 88 
Cale. VII Orbital x 
Carboxyl ion g aS 
Cale. VIII Orbital 2P, | 30, Par of (22) 
Phosphate ion C 1.88 1.60 a 1.88 
Cale. IX Orbital 2p, | 4p, | 2p,, 1-bond fe 
Phosphate ion £ 1.88 | 0.43 1.88 
Cale. X, XI Orbital 2 p,(*)| 4p, (°) 
ADP and ATP C 1.88 0.43 


® For all the O-atoms. 2-Bonds between the P and O atoms. ° For all the P-atoms. 
Table 7. Distances adopted in the calculations (in A). 
The notations for the atoms are given in Figs. 1 and 2. 


Distances in A 
Calculation 


-P 
OL=C C-Og Og-P P-O, | 
I 1.21 1.43 1.60 1.45 1.60 
II 1.21 1.43 1.60 1.55 1.60 
Carboxyl-phosphate IIL 1.21 1.43 1.55 1.55 0.43 
: IV 121 1.43 1.60 1.55 1.72 
Vv 1.23 1.40 1.60 1.50 1:72 
VI 123 1.40 1.60 1.50 0.43 
Carboxyl ion | VII | 1,28 1.28 | 
} 
: VIII i 1.55 1.55 1.60 
mjaephate 19m | IX | 1.85 1.55 | 0.43 
ADP, ATP | 2. A | 1.60% 1.50° | 0.43 


* Distances P-O,, and P-O,,. ° Distances P-O,, P-O, and P-Og. 


boxyl-phosphate 
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but 


Ustru 


B. GRABE, The electron d 


ad UT ag L aA6&GUT nAdU 
eF'0 180= - -18't-|-s90+ LEO | Se0r SET [OS 0+ Po0 

A® 6S — =70% fA9 0S — =9903 ‘a 

eF'0 FLO0-— FLT | SFO+ 2:0 | PL:;0-  FL'I 
e90-— e991 | LZ0+ €2°0 
eF'0 €L0-— ELI | 8FO+  @9°0 | 1€0+ 69°I | 820+ ~ ZL'0 

A® 8T— =20% fae QT — =99Z “0 

eF'0 FLO0— Plt | SFO+ Z9'°0 | FL:0— FLT 
290-— 29° | ee0+ 2L2L°0 
€F'0 SL0= “SLT [°SFO+- 39:0" | SE0+  —S9°L.| Fo:0+ <9L70 

A? 8'1— =9%9¢ the HT — =99S “a 

ar ero | — s90-— sort | 980+ +90 | 890-  89°T 

(6) Jo" O91 | 8SFO- SFT | — _- GSO+ GLO | LLO-  LL'T 
ss0-— SST | 910+ ##8°0 
as ero | — — oso- OST | Se0+ 90 | Ge0+ 89'T | 810+ Z8°0 
(ZZ) Jo a GLI | ¢PO- srt | — _— 840+ G40 | T€0+ 691 | LIO+ £€8°0 
(Go)'30° SLT | OF0— OFT || — — 920+  FL:0 | 930+ FLT | LTO+ ~— €8°0 
ay ero | — — 990-— 991 | 980+  49°0 | 680+ 191 | LIO+ €8°0 
(GG)305-O 2 09° (1P0— TFS | — aa 970+ FLO | Lo0+ ELT | LIO+ £870 
a oor | — — sc0-— sgt | se0+ S90 | €80+ 2L29T | LIO+ &8€8°0 

Ae F1— =99% fA2 OT = = 995 -¥ 
“Fo)ro-qag 45 [4fO)@ “Fol | (“oO (“O)W] (DG (aN | CO) Con] (OO ()N 


uor aynydsoyd pun uor pfiixoguvo ‘aynydsoyd-jixoginp 


49 O~ 
OF 0 — 


€9°0— 
S¢'0—- 


69°0 — 
€€°0 — 


80 — 
FE0— 
FE0— 
660 — 
L160 — 
620 — 
86°0 — 


(70) 0 


*L pure g so[qey, os[e puw suUUUN[OD 04} 4SRT OY} OES ‘sUOTIBTNOTRO YY JO WOT}BIOUINU 
oy} Jo qaodurr oy} ynoqe uoTywUTIOFUT IOYIINJ JOT “GLY pue qqV sed4yz oy} jo soyeydsoyddjod 0, TX pue x pue uor oyeydsoyd oy 03 
XI pue TITA “uor [Axoqreo oy} 03 TTA ‘oyeydsoyd-;Axoqzeo 03 Jeyor [A] SUOTVR[NOTV OY, *Z PUB | “SBI oes ‘suI0ze yy IOF sUOTYRIOU YY IO 


40 - 
OFT 


8ST 
FET 
Lom | 
62° 
LeL 
6e'L 
81 


(70) N 


(ST) puv (LT) 04 Surpaoooe pozr[noyeo ‘, suorjoape-2,, oy} 0} onp (¥)~ pur (¥) Nv *9 29D, 


iA 
TIT 


axXxI=xXI 
ITA 
Tit 


“orep 


216 


ARKIV FOR FYSIK. Bd 15 nr 17 


990- e9°T en OF'0 BRE: T6T |490- LOT [80+ zcr0 |60°0+ TI6T |o9'0+ OF 0 |G9'0—  eo'T 1196 
s90- 9'T |090+ OF 0 |60:0+ T6T |09°0+ OF 0 |e90—- eg'T x 

A® ZZ— =1O% ‘a 
820-— 8'T aes 90 aoe 6L1 |s90-— gor |er0+ L9°0 |IZO+ 6L°T pee €o'0 |sso-— set | IX 
8S0— sot |LF0+ ec°0 |eo0+ set |LFO+ E90 |8e0—  gerT 26 

A® 81— =90% ‘9=_q 
Ig0-— I'L Fie te 990 |600+ TL'T iene 8o'L | cso+ soo [6zo+ Lt |9¢0+4 90 1te0— me Lx 
oso- OST |gs0+ 99:0 |Ts80+ 69T |Ssg0+ 90 |oco-— ogT x 


Ae FIT- =10¢ “vy 


_ SSS eS SE EEE Eee 


?0)® Po)N | (MO (aw |(*0)@ O)N| PO) Po) | (a)O (a) ("0)0 (o)N| (AO (a) | (?0)6 PO)N| “TRO 


217 


B. GRABE, The electron distribution in carboxyl-phosphate 


Table 9. ¢-values (in eV) obtained from z-values according to (10). 


The numeration of the calculations indicates the type of AO’s used and the distances between 
the atoms (see Tables 6 and 7), while the notation A, B ... indicates the 6” values used. The 
Oau values are given in eV. 


Carboxyl-phosphate 


t)) noe Rae ee eee III 
Cale. 
A jorn. |" ohare 
2 066 Al ta(t) 1.4 1.6 “ 
OP 14 =e 1,8.) See 
Orb. 1 14.0 14.0 14.1 14.0) | s>14.0 | —14.0)|" —142) |@—1435 a 1a 
2 13.1 13.5 13.4 13.5 13.3 13.4 13.5 13.7 | —18:8 
3 12.3 12.3 12.3 12.3 12.3 12.3 12.7 12.9 | = 1333 
4 = 3:8 7.9 9.6 7.8 8.1 9.5 9.9 10.1 10.4 
5 6.7 6.3 7.0 6.3 6.4 7p! 7.4 7.5 8.0 
« 
Carboxyl ion Phosphate ion 
VII VIII IX 
Cale a 
AD pis Riad D A A Sie 6) hep 
2056 1.0 = Ld 1.6 2.0 1.0 1.0 1.4 1.6 2.0 
2 bor = - 1.4] — 1.4 1.8 1.8 2.2 
Orb. 1 — 14.0 -| = 14-2 ~ 14.3 14.4 13.6 13:4 | —3;65)5— 138 3) eae 
2 12.3 12.7 12.9 13.3 12.3 12.3 12.7 12.9 13.3 
3 6.9 7.2 7.4 7.8 6.7 =0 9.574" = 99.8. |S OM) ie aes 
Adenosine diphosphate Adenosine triphosphate 
x XI 
Cale. 
Agente boone. ARIS eS D 
it. 1.0 1.4 16| — 20] - 10| — 14] - 16] - 2.0 
2 0op — 14 1.8 1.8 2.2| — 14 1.8 1.8 22 
Orb. 1 ale) || elas) tO 14.0 13.9 13.9 14.1 14.0 
2 13.3 13.5 13.7 13.9 Legion) — 13:6) |= 13:8 ea S-9 
3 = 12.3) — 12.7 12.9 13.3 13.3 13.5 13.74 palo 
4 9.7 9.9 10.1 | —10.3 12.3 12.7 12.9 | —13.3 
5 9.1 9.5 9.7 10.2 T2oel bo Leis ee 20 eel SS 
6 Se Oot 9.9) Ne Oe lel et 
7 = 9.4 |-— 9.7.) = 9.9) = 10:3 
8 55.0) )| — 9.3. |8—=— 054 SO 
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5. Discussion 


The semi-empirical method described above is unfortunately approximate in 
several respects. Firstly, all the necessary empirical data have not been determined 
experimentally so that some of them must be estimated with guidance from other 
existing data. Furthermore there is no quantity that can be determined experiment- 
ally which really corresponds to the 6-integrals. Another point is that it is difficult 
to choose appropriate orbitals for the P-atom and thus for the phosphory] group. 
For this reason, it seemed to be necessary to examine the influence on the results 
of the different parameters in equation (11). 

Two types of variation have been made: (a) Calculations have been performed in 
which the values assumed for the interatomic distances and the AO’s chosen for 
the phosphoryl group were different. In all these calculations (calc. [A-VI A), the 
same values for the 6’s were used, and thus only the non-diagonal elements in (11) 
were varied. (b) The same interatomic distances and AO’s as in IIIA were used but 
the 6-values were varied (calc. IIIB, IITC and III D) and hence all the elements in 
(11) were therefore varied. This latter type of variation was also performed for the 
carboxyl ion, the phosphate ion and the polyphosphates. 

It is seen that the variations of type (a) had only a relatively small influence on 
the results. The characteristics of the distribution of the charges between the atoms 
is the same in all the cases (see Table 8). 

The variations in the 6-values had a greater influence on the charges on the atoms. 
It is seen from Table 8 that the positive charges on the bridge O-atoms in the carboxy]- 
phosphate and the polyphosphates decrease when more negative d-values are chosen 
while the negative charges on the other O-atoms simultaneously increase. In all the 
calculations, however, the charges on the bridge O-atoms are very different from 
the charges on the other O-atoms considered, viz. the other O-atoms in the same 
group as well as those in the simple carboxyl and phosphate ions. These differences 
in charge remain approximately unchanged when the 6-values are varied. For all of 
the d-values used in the present calculations, the bridge O-atoms are positively 
charged and the charges are considerable in all the cases except X D and XID. 

It seems probable that the correct d-values are within the range of values used 
in the calculations. When using the D,-values of Table 4, larger 6-values cannot be 
used in the calculation of y according to (14) since the quantity under the radical 
sign then becomes negative. In other cases, when other D,-values are used, there 
is of course the possibility of using larger 6-values. 

Furthermore, it appears from Table 8 that the charges obtained on C and P 
respectively in carboxyl-phosphate and the polyphosphates are approximately the 
same as the charges on the corresponding atoms in the simple carboxyl and phos- 
phate ions when the same AOQO’s and 6-values are used in the calculations. (Compare 
III, VII, 1X, X and XI for each 6-value.) 

The results of the calculations of the charges due a the *‘m-electrons’’ can be 
summarized as follows. The charges on the C-atom and the P-atoms in the carboxyl- 
phosphate and the polyphosphates are the same as on the corresponding atoms in 
the carboxyl and phosphate ions. The charges on the bridge O-atoms in the carboxy]- 
phosphate and the polyphosphates, however, are very different from the charges 
on the other O-atoms of the same group and from the charges on the O-atoms of the 
carboxyl and phosphate ions, all of the non-bridge charges being much more negative. 
Furthermore it seems very probable that the charges on the bridge O-atoms are 
positive and of considerable magnitude. 
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Fig. 3. An ADP ion attacking carboxyl-phosphate, resulting in the formation of a carboxy]! ion 
and ATP. 


These results seem to give information which is adequate for the purpose of this 
investigation, viz. to gain an insight into the mechanism of the transfer of energy 
in the living cell. It is well known that, when carboxyl-phosphate reacts with other 
molecules, either the P—O bond or the C—O bond can be broken. For example, when 
it reacts with adenosine diphosphate (ADP), the P—O bond is broken and adenosine 
triphosphate (ATP) and carboxyl ion are formed. This is a reaction which requires 
an activated phosphoryl group. In reactions requiring activated acetyl groups, 
however, the C—O bond is broken. 

Using the results of the present calculations, these experimental data can be 
explained in the following manner. Since oxygen is a very electronegative atom, the 
existing distribution of the “z-electrons” undoubtedly causes the c-electrons in the 
o-bonds between the bridge O-atom and the C- and P-atoms to be attracted closer 
to the O-atom than in the simple carboxyl and phosphate ions.! But this is not 
possible without making the C- and P-atoms more positive than in the simple ions. 
Thus, the three atoms in the -C-O—P- group compete for the electrons. This stabilizes 
the molecule and causes the so-called high-energy of the molecule. When the C-atom 
or the P-atom is attacked by a negative ion, the o-electrons of the bond between the 
attacked atom and the bridge O-atom are displaced towards the O-atom and this 
finally results in a splitting of the bond. Let us assume that the P-atom is attacked 
by an ADP ion. This action results in the formation of a carboxyl ion and the phos- 
phorylation of ADP. (See Fig. 3.) Apparently, the attraction from the O-atom on 
the o-electrons of the P—O bond is powerful enough to make the P-atom so positive 
that it becomes able to attract the electrons of ~O—P- in the ADP ion to such a 
degree that this O-atom becomes a bridge atom at the same time as the P—O bond in 
the carboxyl-phosphate is broken and the carboxyl ion is formed. Thus it is the 
powerful attraction from the bridge atom on the o-electrons which makes possible 
the simultaneous formation of a carboxyl ion and a new bridge atom, viz. that in 
ATP. 


1 This view is consistent with results obtained by Fischer Hjalmars [12] in a complete quantum 
mechanical calculation of the electron distribution in the O, molecule. Fischer-Hjalmars found 
that, for the O, molecule, the effect of the z-electron distribution is diminished considerably by 
the o-electron distribution. 
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In connection with the treatment of the synthesis of acylphosphates, Avison [13] 
remarks that acyl-phosphoric acids do not react with orthophosphate ions since their 
negative charge prevents the approach of the attacking phosphate ions. He suggested 
that, in biological systems in similar reactions, the negative charge can be screened 
by chelation or in some other way as a result of the enzymic environment. This 
point of view should be applied to the reactions described in the preceding paragraph. 
In [3, 4], the present author has suggested that there are balanced forces in the car- 
boxyl-phosphate (in the present paper expressed as the competition for the electrons 
by the atoms of the -C-O-P- group) and that a small perturbation of this balance 
by a specific enzyme may render possible a reaction with for example an ADP ion. 
However, according to Avison’s experiments, an attack by a negative ADP ion 
should be sufficient for reaction. Since the carboxyl-phosphate, which is dissociated, 
is protected from the attacks of negative ions by its own negative charge, it is very 
probable that the function of the specific enzyme instead is to render possible the 
attack by screening the charge of the carboxyl-phosphate. In this way the specific 
enzymes may be used to control the nature of the attacking ions and thus the 
anabolism, and it seems probable that the negative charge on the carboxyl-phos- 
phate, protecting against attacks by arbitrary negative ions, makes it especially 
appropriate for energy transfer in the living cell. 

Avison’s explanation of the acylation of orthophosphate ions by carboxylic acid 
anhydrides is different from the present explanation of the phosphorylation of ADP 
ions even though these reactions are similar. He does not take the charge on the 
bridge atom into consideration. As appears from the present investigation, it is very 
probable that this charge is of significant importance for the reaction. 

Recently, Clark et al. [14] have pointed out the necessity of the heterolysis of the 
P-O bond in the phosphorylation of ADP and in similar reactions. 
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APPENDIX 
Table 10. N(j; A) of (16) for each of the lowest space orbitals. 


The notation for the atoms is given in Figs. 1 and 2. The numeration of the calculations indicates 
the types of AO’s used and the interatomic distances (see Tables 6 and 7). 
(i) Carboxyl-phosphate 
Calculations I-VI 


Cale. | Orb. | O« 


Sr? 


Cc 


Op 


Ae Oy1 


A. 269g= —1.0 eV; 2dgp = —1.4 eV 


I 1 0.49 0.41 0.09 0.01 0.00 — 
2 0.05 0.00 0.20 0.32 0.42 we 
3 0.10 0.00 0.54 0.00 0.36 — 
II 1 0.48 0.40 0.10 0.01 -— 0.01 
2 0.05 0.01 0.08 0.35 --- 0.51 
3 0.12 0.00 0.69 0.00 a 0.19 
IIr i 0.44 0.39 0.13 0.03 0.01 — 
2 0.12 0.02 0.21 0.29 0.36 — 
3 0.08 0.00 0.46 0.00 0.46 — 
IV 1 0.48 0.40 0.10 0.01 oo 0.01 
2 0.05 0.01 0.08 0.35 a 0.51 
3 0.12 0.00 0.69 0.00 -— 0.18 
Vv 1 0.44 0.40 0.14 0.02 — 0.01 
2 0.06 0.02 0.08 0.34 — 0.50 
3 0.17 0.00 0.63 0.00 — 0.20 
VI 1 0.40 0.38 0.16 0.03 0.03 — 
2 0.11 0.03 0.07 0.29 0.50 oa 
3 0.16 0.00 0.61 0.00 0.23 oo 
B. 2099= — 1.4 eV; 205,= —1.8 eV 
IIt il 0.51 0.37 0.10 0.01 0.01 
z 0.09 0.01 0.27 0.24 0.39 
3 0.06 0.00 0.47 0.00 0.47 
C. 26o9= — 1.6 eV; 26op= —1.8 eV 
Tit = O54 0.35 0.09 0.01 0.01 
2 0.09 0.01 0.28 0.24 0.38 
3 0.05 0.00 0.47 0.00 0.47 
D. 2699= —2.0 eV; 2pp = —2.2 eV 
TIL 1 0.65 0.32 0.02 0.00 0.00 
2 0.03 0.00 0.38 0.18 0.41 
3 0.02 0.00 0.49 0.00 0.49 
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Table 10 (continued) 


(ii) Carboxyl ion 
Calculation VIT 
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(iii) Phosphate ion 
Calculations VIII and IX 


Orb. | Ox C | Og Calc. Orb. Oz | 1 | Oy; | (O3)y 
A. 2655= —1.0 eV A. 20op= —1.4 eV 
1 0.29 0.42 0.29 VIII 1 0.20 0.38 = 
2 0.50 0.00 0.50 2 68 0.00 = 
B. 2099= —1.4eV A. 20op= = 1.4 eV 
1 0.31 0.38 0.31 IX 1 0.34 0.32 0.34 
2 0.50 0.00 0.50 2 0.50 0.00 0.50 
C. 206¢= —1.6 eV B=C. 2d5p= —1.8 eV 
1 0.32 0.37 0.32 1 0.37 0.26 0.37 
2 0.50 0.00 0.50 2 0.50 0.00 0.50 
D. 205¢= — 2.0 eV D. 2d5p= —2.2 eV 
1 0.34 0.33 0.34 1 0.41 0.19 | 0.41 
2 0.50 0.00 0.50 2 0.50 0.00 0.50 
(iv) Adenosine diphosphate 
Calculation X 
Orb. | O, | ang On | ie | O, 
A. 20op= —1.4 eV 
1 0.1 0.17 0.28 0,17 0.19 
2 0.38 0.15 0.00 0.15 0.35 
3 0.23 0.00 0.57 0.00 0.22 
B=C. 2dpp= —1.8 eV 
1 0.23 0.14 0.25 0.14 0.23 
2 0.38 0.12 0.00 0.12 0.38 
3 0.18 0.00 0.64 0.00 0.18 
D. 2:65, = —2.2:0V. 
1 0.32 0.11 0.15 0.11 0.32 
2 0.41 0.09 0.00 0.09 0.41 
3 0.10- 0.00 0.81 0.00 0.10 
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Table 10 (continued) 


(v) Adenosine triphosphate 
Calculation XI 


Orb. | 0, | P | On; P | 0, | Ors | P | 0, 
A. 26op= -1.4 eV 
1 0.05 0.06 0.17 0.24 0.20 0.17 0.06 0.05 
2 0.25 0.16 0.09 0.00 0.00 0.09 0.16 0.25 
3 0.23 0.10 0.00 0.10 0.23 0.00 0.10 0.23 
4 0.09 0.00 0.23 0.00 0.35 0.23 0.00 0.09 
5 0.14 0.00 0.36 0.00 0.00 0.36 0.00 0.14 
B=C.. 2055 = all SiON 
1 0.06 0.05 0.16 0.20 | 0.25 0.16 0.05 0.06 
2 0.29 0.13 0.08 0.00 0.00 0.08 0.13 0.29 
3 0.25 0.08 0.00 0.08 0.25 0.00 0.08 0.25 
4 0.08 0.00 0.27 0.00 0.31 0.27 0.00 0.08 
5 0.11 0.00 0.39 | 0.00 0.00 0.39 0.00 0.11 
D. 2dop= —2.2 eV 
1 0.10 0.04 0.10 0.15 0.39 0.10 0.04 0.10 
2 0.36 0.10 0.04 0.00 0.00 0.04 0.10 0.36 
io 0.27 0.06 0.00 0.06 0.27 0.00 0.06 0.27 
4 0.04 0.00 0.37 0.00 0.18 0.37 0.00 0.04 
5 0.05 0.00 0.45 0.00 0.00 0.45 0.00 0.05 


Tryckt den 4 april 1959 


Uppsala 1959. Almqvist & Wiksells Boktryckeri AB 
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